Measurements of the Sun's High Latitude Meridional Circulation 
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ABSTRACT 

The meridional circulation at high latitudes is crucial to the build-up and reversal of the 
Sun's polar magnetic fields. Here we characterize the axisymmetric flows by applying a magnetic 
feature cross-correlation procedure to high resolution magnetograms obtained by the Helioseismic 
and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). We focus on 
Carrington Rotations 2096-2107 (April 2010 to March 2011) - the overlap interval between HMI 
and the Michelson Doppler Investigation (MDI). HMI magnetograms averaged over 720 seconds 
are first mapped into heliographic coordinates. Strips from these maps are then cross-correlated 
to determine the distances in latitude and longitude that the magnetic element pattern has 
moved, thus providing meridional flow and differential rotation velocities for each rotation of the 
Sun. Flow velocities were averaged for the overlap interval and compared to results obtained 
from MDI data. This comparison indicates that these HMI images are rotated counter-clockwise 
by 0.075° with respect to the Sun's rotation axis. The profiles indicate that HMI data can be 
used to reliably measure these axisymmetric flow velocities to at least within 5° of the poles. 
Unlike the noisier MDI measurements, no evidence of a meridional flow counter-cell is seen in 
either hemisphere with the HMI measurements: poleward flow continues all the way to the poles. 
Slight North-South asymmetries are observed in the meridional flow. These asymmetries should 
contribute to the observed asymmetries in the polar fields and the timing of their reversals. 

Subject headings: Sun: dynamo, Sun: rotation, Sun: surface magnetism 



1. INTRODUCTION 

The magnetic fields at the surface of the sun 
provide the inner boundary condition for the he- 
liosphere. As such, motions of magnetic flux near 
the solar surface are crucial to studies involving 
the heliospheric magnetic field configuration, ac- 
tive region evolution, the buildup of polar fields, 
and subsequent magnetic reversals. Furthermore, 



surface magnetic flux transport is vital to dynamo 
models used to explain the sunspot cycle itself. 

The axisymmetric flows were characterized us- 
ing magn etic feature tracking for nearly the e ntire 
MDI era (|Hathawav fc Rightmire1l2010l[2^lTh . It 
was found that the meridional flow varied consid- 
erably during that time. Specifically, the merid- 
ional flow speed that led up to Solar Cycle 23/24 
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minimum was much faster than the meridional 
flow speed during the prior minimum. Faster 
meridional flow in the active latitudes (equator- 
ward of ~ 40°) inhibits the cancellation of oppo- 
site polarity magnetic elements across the equator. 
This reduces the imbalance of magnetic polarities 
in the active latitudes which, when transported to 
the poles, leads to weaker polar fields, a weaker 
Solar Cycle 24, and an extended solar minimum 
between Solar Cycles 23 and 24. While this gives 
a credible physical explanation for the peculiari- 
ties of the Solar Cycle 23/24 minimum, it requires 
knowledge of the poleward transport all the way 
to the polar regions. For this reason, as well as 
others, it is important to measure the meridional 
flow to the highest latitudes possible. 

SDO was lau nched in February, 2 010. On 



board was HMI ( Sche rrer et al. 112012ft . HMI is 



a follow-on and more capable version of MDI. 
HMI magnetograms, with a size of 4096 2 , have 
four times the spatial resolution of MDI full-disk 
magnetograms. HMI magnetograms are contin- 
uously available with a cadence of 45 seconds, 
rather than the 96 minutes of MDI. HMI mag- 
netograms averaged over 720 s econds are virtu - 



ally unaffected by the p-modes (|Liu et al. 1 12012) 



These advances make HMI data ideal for the con- 
tinuation of the correla t ion trackin g analysis of 
Hathaway fc Rightmirel (|2010l l201lh and ideal for 
extending the measurements to higher latitudes. 

The meridional flow at high latitudes has sig- 
nificant consequences for models of the Sun's mag- 
netic dynamo. Models of the magnetic flux trans- 
port at the Sun's surface (Ivan Ballcgooi ien et al 



19981 : ISchriiver fc Title ifcoOll : IWang eta!. I l2009h 
have employed meridional flow profiles which ei- 
ther stop completely at 75° latitude or quickly fall 
to zero before entering the polar regions. Flux 
Transport Dynamo models depend critically on 
the strength and structure of the meridional circu- 
lation and, with assumptions about the meridional 
circulation, have been used to pre dict the ampli 



tude and timing of Solar Cycl e 24 (|Dikpati et al. 
20061 : IChoudhuri et al. 112007ft . The presence or 



absence of counter-cells in the meridional flow at 
high latitudes has been shown to alt er the length of 



the s unspot cycle in these models (jDikpati et al. 
2010ft . 



ited to lower latitudes. Komm et al. ( 1993ft 
limited thei r measur e ments to latitudes less 
than 52.5°. iMeunierl dl999ft included measure- 
ment to 70° while lHathawav fe Rightmire I (|2010h 
stopped at 75°. Measurements of the merid- 
ional flow using the methods of helioseismol- 
ogy have been limited to latitudes below 50° 
dGonzalez Hernandez et al~ 2010t Basu fc Antia 
l201f)ft . While direct Doppler measurements can 
conceivabl y measure th e meridional flow right to 
the poles ([Ulrich 112010ft . these measurements are 
subject to systematic errors from the Convective 
Blue Shift signal (an apparent blue shift of spec- 
tral lines due to the correlation between emergent 
intensity and radial velocity in convective flows 
at the surface). This paper seeks to measure the 
axisymmetric flow velocities, especially the merid- 
ional flow, at the highest latitudes possible using 
magnetic feature tracking with HMI data. 

2. ANALYSIS 

The HMI 720 second magnetogram data were 
selected for this investigation due to their re- 
duction in nois e, particularly near the limb 
(|Liu et al. 1 120121 ) . While spatial resolution is im- 
paired at higher latitudes, this only results in 
noisier measurements at the higher latitudes (an 
effect that is mitigated by the high resolution of 
HMI) without inducing a net flow signal. For this 
study we chose HMI 720 second full-disk line-of- 
sight magnetograms obtained every hour on the 
hour. We account for the radial nature of these 
magnetic elements by dividing the line-of-sight 
signal by the cosine of th e heliographic angle from 
disc center as done by lHathawav fc Rightmire 



(l2010Ll201lft . The resulting hourly magnetograms 
were then mapped into heliographic coordinates, 
while excluding data within three pixels from the 
limb. Corrections were made for the ephemeris 
error in the orientation of the Suns rotation 



Previous measurements of the meridional mo- 
tion of the magnetic elements have been lim- 



axis (iBeck fc Giles Il2005t lHathawav fc Rightmire 
2010h . Figure 1 shows an example of a mapped 
image from a 720 second magnetogram. 

Measurements of the axisymmetric flows re- 
quire averaging over all longitudes during a full 
rotation of the Sun. Here Carrington Rotations 
2096-2107 were chosen from the ~ 1 year of over- 
lap between the HMI and MDI observations. Par- 
ticular attention was paid to the systematic and 
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Fig. 1. — A 720 second averaged magnetogram 
from 9/2/2010 at 10:00:00 U.T. mapped into heli- 
ographic coordinates from pole to pole and ±90° 
from the central meridian and with a resolution 
of 4096 2 . For this date and time the North po- 
lar region is well sampled while the South is not. 
An example of a 41x2401 pixel strip used in the 
cross-correlation analysis is marked for reference. 

statistical errors and the consequent latitudinal 
limits for the two data sources. The measurement 
method is based on geometry. Errors in the size, 
shape, or orientation of the images can lead to sys- 
tematic errors in the desired flows. Using simulta- 
neous measurements from these two instruments 
allows us to identify these errors. 

The 4096x4096 magnetic maps were divided 
into 41x2401 pixel strips, one for every 10 th pixel 
latitude position. (One of these strips is shown in 
Figure 1 . ) The HMI strips correspond to the same 
area as the 11x601 pixel MDI strips. In both cases, 
the latitudinal resolution of the cross-correlations 
is limited by the strip width to about one half 
of a degree. Each strip was then cross-correlated 
with corresponding strips from 8 hours later (no 
measurements were taken if the shifted strips ex- 
tend ed beyond the mapped data). Previous test- 
ing (jHathawav fc Rightmire 11201 If) indicates that 
this cross-correlation method does in fact measure 
the axisymmetric motions of the magnetic element 



pattern. The diffusive effects of the random mo- 
tions of the magnetic elements due to the nonax- 
isymmetric convective flows have little or no sys- 
tematic effects on the measurement results. Mea- 
surements were obtained for ~ 600 image pairs 
during each 27-day rotation of the Sun. These 
measurements were then averaged over the entire 
rotation and the RMS variations about those av- 
erages were calculated. 

The cross-correlation analysis was further im- 
proved by incorporating a Forward-Backward 
technique. In the Forward step (previously used 
exclusively) a strip centered on the central merid- 
ian is taken from the initial map and cross- 
correlated with shifted strips from the later map. 
In this case, the displacement is typically prograde 
and poleward. In the Backward step, a strip cen- 
tered on the central meridian is taken from the 
later map and cross-correlated with shifted strips 
from the initial map. In this case, the displace- 
ment is typically retrograde and equatorward. 
The correlations were then combined and the 
peaks in the subsequent cross-correlations were fit 
to parabolas with the position of that peak used 
to determine flow velocities in both latitude and 
longitude. By applying this Forward-Backward 
technique, we are able to double our statistics 
to further reduce the noise. Moreover, any sys- 
tematic effect that might mimic flow toward or 
away from disk center is then canceled in the 
meridional flow measurements by this Forward- 
Backward technique. 

Measurements were attempted at latitudes up 
to 85° in both the North and South but were lim- 
ited to latitudes less than ~75° when the tilt of 
the Sun's rotation axis was unfavorable. Analysis 
equatorward of 30° requires masking out active 
regions, whose flows are not representative of the 
surface axisymmetric flows. Pixels with B > 1000 
G and the adjacent pixels were masked out in or- 
der to exclude flows representative of active re- 
gions. 

3. RESULTS 

Axisymmetric flow velocities from each Car- 
rington Rotation were averaged over the Carring- 
ton Rotations investigated in this study (2096- 
2107) using the inverse of the standard errors 
squared as weights. The HMI average flow veloc- 
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ity profiles were then plotted on top of the average 
flow velocities obtained using MDI data from the 
same time periods. Results obtained with HMI 
have significantly less noise than those from MDI, 
especially near the poles. 

The HMI meridional flow results were found to 
be ~ 2 — 3 m s _1 more northward than the MDI 
results. This offset can be explained by a rota- 
tion of the HMI imagi ng system with re spect to 



MDI, as suggested by iLiu et al~l (|2012t) . (Note 



that the MDI images were corrected for a 0.19° 
rotation relative to the Sun's rotation axis.) Ap- 
plying a least squares fit to the difference between 
the two profiles indicates a counter-clockwise ro- 
tation of HMI by 0.075° ( or a velocity correction 
of ^2.5cos(£?) m s -1 , where B is the heliographic 
latitude). The corrected HMI average flow veloc- 
ity profiles are plotted along with the MDI profiles. 
The differential rotation is shown in Figure 2 and 
the meridional flow is shown Figure 3. 
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Fig. 2. — The Differential Rotation from Carring- 
ton Rotations 2096-2007. The MDI result is plot- 
ted in black and the HMI result is plotted in red. 
The five sigma error ranges are indicated by the 
dotted lines. 

The MDI and HMI differential rotation equa- 
torward of 55° are in very good agreement. There 
are however, small systematic differences (~ 1 — 2 
m s _1 ), with larger velocities in HMI. Poleward of 
55° HMI is ~ 10 — 20 m s _1 more retrograde rela- 
tive to the rotating frame of reference than MDI. 
This may be attributed to a reported elliptical dis- 
tortion of the MDI images ( Korzennik et al. Il993h 
which was not completely accounted for in our 
mapping due to the fact that it is not well char- 
acterized. The statistical errors in the HMI data 
(dotted red lines in Figure 2) indicate that precise 
measurements can be made to ±85° latitude. 
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Fig. 3. — The Meridional Flow from Carrington 
Rotations 2096-2007. a) The corrected HMI re- 
sults are shown in red. Results from MDI are 
shown in black. The two sigma error ranges are 
indicated by the dotted lines, b) The poleward 
velocities in the North (red) and South (blue) are 
plotted to highlight the North- South asymmetry. 



HMI measurements of the meridional flow also 
agree well with those obtained at latitudes up to 
~ 55° (Figure 3a), but are systematically slower 
by ~ 1 — 2 m s _1 . Again, this can be attributed 
to a possible elliptical distortion of MDI images. 
At higher latitudes HMI tells a different story 
than MDI. The MDI measurements suggested a 
counter-cell in the North (equatorward flow above 
~ 60°) but flow to the pole in the South. The HMI 
measurements have no indication of counter-cells 
in either hemisphere for this time interval. The re- 
sults do show a slight, but potentially important, 
North-South asymmetry. A faster poleward flow 
is seen in the South through the active latitudes 
(from the equator to ~ 40°) and in the North at 
higher latitudes. The errors in the HMI meridional 
flow data (dotted red lines in Figure 3a) indicate 
that many of fluctuations in the meridional flow 
profile are actual features, rather than noise. The 
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noise level is low enough that the meridional flow 
is well determined at 85° using HMI data. 

4. CONCLUSIONS 

We find that the improved data from HMI 
can extend our measurements of the axisymmct- 
ric flows, differential rotation and meridional flow, 
to much higher latitudes - 85° or more. During 
this analysis, measurements were restricted to 85°. 
While noise levels are higher at the poles, it fea- 
sible to expect that measurements could be ob- 
tained all the way to the pole, particularly for 
Carrington Rotations in which the sun is tilted 
7° toward or away from the earth. Our measure- 
ments of these flows during the ~ 1 year of overlap 
between the MDI and HMI instrument operations 
clearly show that the poleward meridional flow ex- 
tends to the poles. We do not find evidence for any 
polar counter-cells as was indicated in the highest 
latitude measurements from MDI. 

Recently IZhao et al. ( 2012 ) found and cor- 



rected a systematic error in the meridional flow 
measurements made with time-distance helioseis- 
mology. Their comparison of their corrected 
meridional flow profile with the contemporaneous 
meridional flow measurement s from magnetic ele- 
ment motions seen with MDI (jHathawav &: Rightmire 
2010) a nd fr om direct Doppler measurements 



show good agreement to latitudes 
of 50 — 60° but not much agreement at higher 
latitudes. The high latitude measurements we 
show here are largely in agreement with the he- 
lioseismology results but further out of line with 
the direct Doppler measurements (which show 
counter-cells above 60° in each hemisphere). 

We also find slight, but potentially important, 
North-South asymmetries in the meridional flow 
profile. The poleward flow is faster in the South in 
the active latitudes and faster in the North in the 
polar latitudes. Both of these asymmetries may 
help to explain the observe d North- South asym - 

" (120121) 



Shiota et al. 



metry in the polar fields, 
observed a faster decline in magnetic flux of the 
North pole between 2008 and 2012 than was seen 
in the South. The trend in the North suggests 
an imminent polar field reversal. While some of 
this trend may be due to the fact the the northern 
hemisphere was more active during this period, 
the meridional flow asymmetry we observe should 



also contribute. 

The Sun's polar field reversals are produced by 
the poleward transport of opposite polarity mag- 
netic flux from the active latitudes (opposite to 
that of the pole at the start of the sunspot cy- 
cle). While active regions (sunspot groups) have 
a balance of both magnetic polarities, in each 
hemisphere the opposite polarity is systematically 
found at higher latitudes (Joy's Law). Some of 
the lower latitude, like polarity, magnetic elements 
can cross the equator to cancel with similar (but 
reverse) polarity elements from the other hemi- 
sphere and leave behind an excess of the higher 
latitude, opposite polarity, magnetic elements. A 
fast meridional flow through the active latitudes 
inhibits this cross-equator cancellation and leaves 
a smaller excess of the opposite polarity for sub- 
sequent transport to the poles. Thus, the faster 
meridional flow in the southern active latitudes 
should result in a slower erosion of the South po- 
lar field. In addition, the fast poleward flow at 
high latitudes in the North should accelerate the 
erosion of the North polar field. 

This poleward flow above 75° latitude is, 
nonetheless, problematic for most models of the 
magnetic flux transport in the near surface layers 
(Ivan Ballegooiien et al. 19981: Schriiver fc Title 



200 li I Wang et al. I I2009T ). In these models this 



high latitude poleward flow tends to produce po- 
lar magnetic fields that are too strong and too 
highly concentrated at the poles themselves. Our 
observations may require further adjustments to 
those models. 

The ultimate consequences of this meridional 
flow and its asymmetries will require continued 
measurements of the (variable) meridional flow to 
high latitudes along with improved models of the 
flux transport. 
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